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Recently, the large CP asymmetries in B — ¥ n ± TT + n~ decays were found by the LHCb Collab- 
oration to localize in the region m^ +7r _ < 0.4 GeV 2 . We find such large localized CP asymmetries 
may be due to the interference between a light scalar and p°(770) intermediate resonances. Con- 
sequently, we argue that the distribution of CP asymmetries in the Dalitz plots of three-body B 
decays could be very helpful for identifying the presence of the scalar resonance. 
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Recently, the LHCb Collaboration found clear evi- 
dence for direct CP violation in some three-body de- 



cay channels of B mesons such as B d 



"7T and 



B ± tt+tt'K^ [l|, 0. Intriguingly, large direct CP 
asymmetries wrere found in some localized phase spaces 
of the two decay channels. For B ± — > -k+tt^-k^ 1 , the 



YM3 —> M1M2M3, respectively, and they take the form 

M x = —(sis - §i 3 )(T x + P X e iSx e^), (3) 
sx 

My = —{T Y + Pye^e^). (4) 
s Y 



CP asymmetry in the region 



9 

7T~ 7 



high 



> 15 GeV z 



low < ^-4 GeV 2 and l n the above two equations, Sij = 1,2,3) is the in- 
variant mass squared of mesons M, and Mj, gx(Y) is the 
effective coupling constant for the strong decay X(Y) —> 
M1M2, s X (y) is the reciprocal of the propagator of X (Y) 



A 



CP 



-0.622 ± 0.075 ± 0.032 ± 0.007, 



(1) 



while in the region 
2 

m rr+rr- high 

observed [J] . 



! + _ , < 0.4 GeV^ and 

T+7T low 

< 15 GeV 2 , no large CP asymmetry was 



In this paper, we will show that the localized large 
CP asymmetry may arise from the interference between 
intermediate p° and another scalar meson nearby in the 
three-body decays. 

It is known that the scalar resonance is very difficult 
to identify because of its large width. In the following, 
we will show that the localized CP asymmetries could 
be very helpful for identifying a scalar resonance which 
interferes with the vector one nearby. We will consider a 
B meson weak decay process, B — > M1M2M3, where Mj 
(i = 1,2,3) is a light pseudoscalar meson. If this process 
is dominated by a resonance X in a certain region of its 
Dalitz plot, then it will be very difficult to tell whether 
another resonance exists close to X. We assume that X 
is a vector meson, the possible resonance Y nearby is a 
scalar meson, and both X and Y decay to M1M2. The 
amplitude for B M1M2M3 around the Y resonance 
region can be expressed as 



M = Mx + M Y e 1 ' 



(2) 



which takes the form S12 - 



T 



' m X(Y) +^ 



X(Y) 



/ Sl2 r X(Y)(si2) II, 

and Px(Y) are the tree and the penguin ampli- 
tudes for the decay B -> X(Y)M 3 , f x = T x /(e* ■ p B ) 
and Px = Px/{z* ' Pb) with e being the polarization 
vector of the meson X [6j, Sx,y are the relative strong 
phases between the tree and the penguin amplitudes, <j> 
is the weak phase, and §13 is the midpoint of the allowed 
range of si 3 , i.e., §13 — (si 3imax + si 3im i n )/2, with s 

13, max 

and Si3 jln i n being the maximum and minimum values of 
S13 for fixed S\2- One can check that 



1 



m M- s 12 



(m 2 Mi -m 2 M2 Xm%-m 2 M3 ) 



S12 



, (5) 

The second term in Eq. ([5]) is small compared with the 



first one, since usually (m 



2 

Mi 



) <C m x . Throughout 



this paper, we will denote the momentum, the mass, and 
the decay width of a particle X by px, mx, and Tx, 
respectively. 

As aforementioned, we will focus on the region around 
the Y resonance, i.e., my — Ai < y^sli < my + A2, 
where Ai and A2 are of the order of Ty. We also require 
that my — Ai > mx + Tx (if my > mx) or my + A2 < 
mx — Tx (if my < mx), so that these two resonances 
have competitive contributions in this region. 

For the region of phase space (denoted by lo) where 
the two amplitudes Mx and My are competitive, the 
direct CP violation parameter is found to be 



where 6 is a relative strong phase, Mx and My are 
the amplitudes for B -» XM 3 MiM 2 M 3 and B 



where 



(6) 
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7 CJ 



T| + Ty + P| + + 2 cos fepx cos S x + T Y P Y cos 5 Y 



>12U*13 



T X P Y sin(<5y + 5) + T Y P X sin(<5 x - 5) 



(7) 

(8) 
(9) 



A%=2f As 12 As 13 [f x f Y cos8 + P X P Y cos(<5 x - 8 Y - 6) + cos<j>[f x P Y cos(8 Y +6)+ f Y P x coa(5 - S x )]}, (10) 



with 5 = 5 + arg(sx) — arg(sy), and 
9x 



T 



x 



\*x\ 

9Y 

\*y\ 



r (si3 - si 3 )Tx, 



Jy — i — t-ly, 



(11) 
(12) 



and similar definitions for Px and Py. 

From Eqs. ^ and (U]), one can easily check the fol- 
lowing relations, 



M x {s 13 ) = -M x (s 13 ), 
M Y (s 13 ) = M Y (si 3 ), 



(13) 
(14) 



where s\ 3 — 2s\ 3 ~s\ 3 . These relations allow us to divide 
naturally the region around the Y resonance into two 
parts: f2 and Q, where fl is for s± 3 > §i 3 and O is for 
S13 < S13. From Eqs. (fT3")) and (|14l) . we can derive the 
following relations between Q and Cl phase spaces: 



Besides the CP violation in Eq. 
quantities 



(15) 
(16) 

we define four other 



R± = 



w+ = 



(N n ± N n ) - {N n ± N n ) 
(N n ± N Q ) + (7V fi ± N n ) ' 
(jV n - iV n ) ± (7V fi - A> n ) 
(7V^ + iv n ) ± (7V° + N n ) ' 



(17) 



(18) 



where all the TV's (TV's) are the event numbers of B — >• 
M X M 2 M 3 (B -> M1M2I3) in the corresponding phase 
space. With Eqs. 

(19) 
(20) 
(21) 

Note that R- is independent of the weak phase <j> and 
|i2+| < 1 by definition. So far, we have six quantities: 
Aq P , A^p, R±, and W±, but only three of them are 



3J) and 


US]), one 


P± = 




W+ = 




W- = 





r 



independent. Alternatively, the CP violations in phase 
spaces f2 and O read 

(22) 



1 + R+ 
1-R- 



1 - p 



(23) 



+ 



'CP' 



One can see that the CP asymmetries in these two 
regions can be very different because of the existence of 
the antisymmetric terms A± under the interchange of 
f2 and f2. These antisymmetric terms exist because the 
two resonances X and Y have different spins. If both X 
and Y have the same spin, then A± = A^_ , and one would 
observe that the CP asymmetries in the two regions equal 
each other. One may argue that we cannot exclude the 
possibility that the CP violations may be the same in 
phase spaces f2 and Cl even if X and V are vector and 
scalar mesons, respectively. This is indeed true, and the 
CP asymmetry difference between f2 and cannot be 
used as a probe of the scalar resonance in this situation. 
However, both P_ and R + become good probes. The 
nonzero values of P_ and i?+ will imply the presence of 
the scalar resonance Y. One can check that if Y is a 
vector resonance, then both R + and i?_ equal zero. 

Furthermore, there is also an alternative criteria that 
can be used to identify the resonance of Y. Since the 
amplitude M x becomes very small when s\ 3 is close to 
S13, the amplitude M. Y will be dominant over M Xl and 
then one should observe a larger density of events when 

512 ~ m Y than when S12 ~ m x , on the condition that 

513 is close to S13 0]. 

We have used the transverse approximation for the 
propagator of the vector meson X. The numerator of 
the propagator of X is <? M „ — k^k v /si2 (up to a phase 
factor) with k = pM t + Pm 2 - This has a different off- 
shell behavior from the propagator for a pointlike vector 
particle, — k^ku/m^. In fact, since hadrons are not 
pointlike particles, one inevitably confronts this kind of 
ambiguity when dealing with vector mesons. If we in- 
stead use the latter form of the propagator for the vector 
resonance, we should add to §i 3 in Eq. (jll[) with a term 
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When s 13 is far away from s 13 , this term is small com- 
pared with S13 — S13. It only becomes comparable with 



3 



sis ^ S13 when S13 is close to §i 3 . However, in this case, 
M.x is small compared with .My. Therefore, we are free 
to use the transverse approximation for the propagator 
of the vector meson. 

We want to mention the following two special cases: 

Case 1 : Both 5x and Sy are very small, but S is not 
small. In this situation, both and are small and 
can be neglected safely. One would observe that Aq P 
and Ap P have opposite signs. 

Case 2: All the three strong phases 5x, <5y, and 
S are very small. In this situation, the CP violation 
parameters in both regions will be very close to zero. 
Then, one cannot identify the presence of Y just through 
the measurement of CP violation parameters. However, 
one can still identify the presence of Y by measuring 
R + . The nonzero value of i?+ indicates the existence of 
B ->■ YM 3 M 1 M 2 M 3 . 

In the above discussion, we have assumed that X and 
Y are vector and scalar mesons, respectively. One can 
arrive at a similar conclusion by reversing their spins. 
Our analysis can also be generalized to situations when 
both X and Y have arbitrary spins. If X is a resonance 
with spin J, the corresponding amplitude Aix would 
be proportional to (si 3 - s^)(s 13 - sf 3 ') ■ ■■ (s 13 - sfo), 

where s^ 3 , s^ 3 , . . . , s^iP ^ e within the allowed range of 
S13. Take X as a tensor meson (Y still a scalar meson) for 

example. In this situation, Aix oc (s\ 3 — S13 )(si3 — s^ 3 ), 
where s^ 3 = S13 — A^/v^ and s^ 3 — §13 + A\ 3 /V3 with 
A13 = (si 3!inax - Si3,min)/ 2 - One would observe that 
there is a large difference of CP asymmetries between 
the middle part (s^ < s\ 3 < §i 3 ) and the other two 

parts (si 3 < and S13 > s^ 3 ). 

Now we are ready to show that the large localized CP 
asymmetries observed by LHCb in B ± — > 7r ± 7r + 7r _ can 
be interpreted as the interference of p° and /q(500). The 
LHCb Collaboration found that for B ± — > 7r ± 7r + 7r~, the 
dominant resonance is the vector meson p°(770) In 
the region sl < 0.4 GeV 2 , there is a large difference of 
CP asymmetries between the upper (sh > 15 GeV 2 ) and 
the lower (sh < 15 GeV 2 ) parts. In the following, we will 
denote these two parts by W and fi', respectively. Note 
that 15 GeV 2 is very close to sh, which is about 14 GeV 2 
for < 0.4 GeV 2 . According to the above analysis, one 



immediately concludes that there is a resonance with spin 
lying in the region s L < 0.4 GeV 2 . From PDG [|[, we 
know that this particle could be /o(500). In the following, 
we will show that by including the amplitudes for B ± — > 
70(500)^* n + TT--K ± and B ± p°n ± tt+tt-tt*, 
the observed CP violation behavior can be understood. 

We assume that the two amplitudes of B ± — > 
70(500)^ ->■ Tr+Tr-TT* and B ± -> p ^ -> tt+tt-tt* 
are dominant for sl < 0.4 GeV 2 . They can be expressed 




FIG. 1. Allowed region for d\ and 5\. If plotted in polar 
coordinate system, one can find the allowed region is actually 
a circular ring crossing the origin. 

as 

M = inZLfa s h ) Mb ;^-^ (24) 
s p o e* ■ p B 

M fo = ^M B -^- fo , (25) 

where /o represents /o(500) and Mb-^-r- p an d 
M.B--*ir- f are the amplitudes for B~ — > ir~ p and 
B~ — >• 7r _ /oi respectively. 

With the effective Hamiltonian for the weak transition 
b — > qqd 0], one can obtain the decay amplitudes for 
B — > pir and B — > /o(500)7r, which can be expressed as 
(a common factor Gp/y/2 has been neglected) 



M, 



- p0 = V ub V: d [a 2 X^P°^ + ai Xi B ^-^} - V tb V t 



Id 



04 + aio + 



—die lSl - 2\ — + a&)m ^ 



(m d + m u ){m b + m u ) 



-a 4 + -a 7 + -a 9 + -a w 



x (B-p°,n-) 



X (B-n-,p°) 



(26) 



M B -^ fo „- = V ub V* d a 2 X^ B fo ' w ^ ~ V tb V t 



Id. 



a 4 + aio 



TUB 



doe 1 



(a 6 + a 8 )m 2 



(m u + m d )(m b + m u ) 



X (B-f ,7T-) 



(27) 



where all the Oi's are built up from the effective Wil- son coefficients C^s, and take the form ai = C[ + 



4 



C' i+1 /N c for odd i and a, = C{ + C' l _ 1 /N c for even 
i, the notation X for matrix elements is borrowed 
from Ref. [lfj- For example, P ° : ' T ) is defined 

as (ir~\(du)v-A\0}(p°\(ub)v-A\B~). These matrix ele- 
ments can be parametrized as the products of decay con- 
stants and form factors. For numerical results, we use 
Ff-^(0) = 0.25 and A^ p {0) = 0.28 We also sim- 
ply set F B ^ fo (Q) = 0.3. 

Terms containing d\ and d 2 come from annihila- 
tion terms, which are proportional to X^ B - p * ) or 
X( B '/"' r ). Usually, annihilation terms are suppressed 
by at least a factor ^QCD/^bj so that one can ne- 
glect them safely. However, there are also annihilation 
terms that are enhanced by a chiral factor, m 2 B /[(rrib ± 
m u )(m,d ± m u )]. This kind of term should be taken into 
account with proper parametrization. According to our 
parametrization, di and d 2 should be, at most, order 
one. Because of multiple soft scattering, annihilation di- 
agrams may also give rise to strong phases. This explains 
the appearance of Si and 8 2 ■ For the effective Wilson co- 
efficients, we will adopt the set of coefficients in Ref. fl^ . 

We have the following five free parameters: 8, di, Si, 
d 2 , and S 2 . Since d\ and 81 are related to the chiral en- 
hancement, this makes them potentially sensitive to the 
branching ratio of B + — > p°Tr + . Thus, these two parame- 
ters can be constrained by the experimental data for the 
branching ratio of B + —> p°ir + . We use the following 
experimental data to determine the allowed region for di 
and 81 y|: 

BK(B+ p°tt+) = (8.3 ±1.2) x 10~ 6 . (28) 

The results are shown in FIG. [TJ 

For given allowed values of di and Si, we should deter- 
mine the allowed regions for the other three parameters 
with the aid of the data, 

A% P = ±0.62 ±0.10, (29) 
a cp = -0.05 ± 0.05. (30) 

In Tablefl] we show the allowed regions of 8, d 2 , and 8 2 for 
given values of d\ and Si. Note that the allowed regions 
of these three parameters are in fact correlated. What 
we show in the table are actually the largest ranges. The 
correlated allowed region of these parameters is a subset 
of the direct combined region shown in Table |U 

The change of input parameters may change the al- 
lowed regions of the parameters shown in Table HI but 
it does not change the conclusion that the large CP 
asymmetry difference between phase spaces fl' and tt' 
is caused by the interference of p° and Jo (500). We also 
anticipate that R± should be nonzero, and this can be 



checked by the data very easily. Because Aq P (so that 
Agp) is very small, we also predict that i?_ is a little bit 
larger than 1. 

We confronted two resonances during our calculations, 
p°(770) and /o(500). The masses and total decay widths 
of these two resonances in our numerical calculation are 
(in GeV) [H 

™V(77o) = 0.775, r p o (770) = 0.149, 

«Vo(5oo) = 0.500, r /o(500) = 0.500. 

Since the nature of /o(500) is not known yet, our numeri- 
cal calculation here should be regarded as an estimation. 
We also used the factorization hypothesis during our cal- 
culations of amplitudes corresponding to the two inter- 
mediate resonances, p° and /o(500). Since the p° meson 
is not on the mass shell, the calculation with this hypoth- 
esis is clearly not accurate. However, since the interested 
area of the phase space is not far away form the p° mass 
shell, the factorization hypothesis is still good enough for 
an estimation. 

In summary, we have shown that the interference of 
two intermediate resonances with different spins can re- 
sult in a CP violation difference in the corresponding 
phase space, which can be used as a method to identify 
the scalar resonance that is close to a vector one. With 
this method, we show that the recently observed large 
CP asymmetry difference in — > ■k ± -k + tt~ decays lo- 
calized in the region m„+^- < 0.4GeV 2 indicates the 
existence of a scalar resonance, which can be identified 
as / (500). 



TABLE I. Allowed regions of S, 82, and t?2 with given values 
of di and Si. 



(di,5i) 


5 


82 


d 2 


(0.7, 260°) 


(3°, 178°) 


(-2°, 36°) U (114°, 153°) 


(0.2,0.6) 


(0.2,190°) 


(-31°, 54°) 


(-1°,26°) 


(0.1,0.4) 


(0.2,330°) 


(141°, 209°) 


(133°, 153°) 


(0.1,0.4) 
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gz 
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